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Abstract. In this paper we provide characterizing properties of TDI systems, among 
others the following: a system of linear inequalities is TDI if and only if its coefficient 
vectors form a Hilbert basis, and there exists a test-set for the system's dual integer pro- 
grams where all test vectors have positive entries equal to 1. Reformulations of this provide 
relations between computational algebra and integer programming and they contain Ap- 
plegate, Cook and McCormick's sufficient condition for the TDI property and Sturmfels' 
theorem relating toric initial ideals generated by square-free monomials to unimodular 
triangulations. We also study the theoretical and practical efficiency and limits of the 
characterizations of the TDI property presented here. 

In the particular case of set packing polyhedra our results correspond to endowing the 
weak perfect graph theorem with an additional, computationally interesting, geometric 
feature: the normal fan of the stable set polytope of a perfect graph can be refined into a 
regular triangulation consisting only of unimodular cones. 



1. Introduction 

A restricted draft concerning an earlier stage of this research has been reported in the 
IPCO 2007 conference proceedings |28j . 

Let A = [ai a2 • • ■ a n ] £ Z rfxn and assume that A has rank d. With an abuse of notation 
the ordered set of vectors consisting of the columns of A will also be denoted by A. For every 
a C [re] := {1, . . . , re} we have the d x \o~\ matrix A a given by the columns of A indexed by 
a. Let cone(^4), ZA and ~NA denote the non-negative real, integer and non-negative integer 
span of A respectively and assume that "LA = 7L d . 

Fixing c £ IR n , for each b € W 1 the linear program (or primal program) LP J 4 )C (b) and its 
dual program DPA, c (b) are defined by 

LP^ jC (b) := minimize {c • x : Ax = b, x > } 

and DP^ c (b) := maximize {y • b : yA < c}. Let Pb and Q c denote the feasible regions 
of LP J 4 jC (b) and DP J 4 jC (b) respectively. Note that the linear program LP^ ;C (b) is feasible if 
and only if b 6 cone(^4). We refer to Schrijver |33j for basic terminology, facts and notations 
about linear programming. 

The corresponding integer program is defined as 

IPA, c (b) := minimize {c • x : Ax = h, x <G N n }. 
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We say that c G W 1 is generic for A if the integer program IP^^b) has a unique optimal 
solution for all b G NA If c is generic then Q c ^ and each linear program LP J 4 c (b) also 
has a unique optimal solution for all b G cone(^4) but the converse is not true in general. 
However, the converse clearly holds for TDI systems: the system yA < c is totally dual 
integral (TDI) if LP,4 jC (b) has an integer optimal solution x G N n for each b G cone(A)nZ rf . 
In other words, the system yA < c defining Q c ^ is TDI exactly if the optimal values 
of LP A , c (b) and of IP^cO) coincide for all b G cone(A) n Z d . (If Q c / 0, LP and DP are 
both feasible and bounded.) This is a slight twist of notation when compared to habits in 
combinatorial optimization: we defined the TDI property for the dual problem, in order to 
be in accordance with notations in computational algebra. 

Totally dual integral (TDI) systems of linear inequalities play a central role in com- 
binatorial optimization. The recognition of TDI systems has been recently proved to be 
coNP-complete by Ding, Feng and Zang [3] and this result has been sharpened to the recog- 
nition of explicitly given systems with only — 1 coefficient vectors, and where the defined 
polyhedron has exactly one vertex (Hilbert basis testing by Pap [30]). Graph theory re- 
sults of Chudnovsky, Cornuejols, Liu, Seymour and Vuskovic [9] allow one to recognize TDI 
systems with — 1 coefficient matrices and right hand sides. However, solving the corre- 
sponding dual pair of integer linear programs (including the coloration of perfect graphs) 
in polynomial time with combinatorial algorithms remains open even in this special case. 
The fixed dimension case has been solved long ago [7], whereas the fixed codimension case 
only recently [13] , [22] , the first using a generalization of integer programming, the second 
computer algebra, both starting from a characterization of Hilbert bases in |35j . 

A particular case where the recognition of TDI systems is still open, is the case of generic 
systems (Problem \3.2\i . which is slightly more general than the perfectness test (detecting 
perfection in a graph), and could be a possible start for an alternative, simpler algorithm 
for the latter. This paper wishes to contribute to testing TDI in cases that occur in integer 
programming and combinatorial optimization, and which usually do not belong to the 
extremities that have been understood so far. 

In Section [21 characterizing properties of TDI systems are provided. Some of these 
properties involve tools from combinatorial optimization, some others from computational 
algebra. Section [3] specializes these results to integral set packing polytopes. Finally, Sec- 
tion [3] exhibits the possible alternatives and their relative efficiency for recognizing TDI 
systems. The remainder of this introduction is devoted to explaining the main results and 
providing some of the necessary background. 

A collection of subsets {a±, . . . , at} of [n] will be called a regular subdivision of A if there 
exists c G W 1 , and zi, . . . ,Zj G such that z.; • &j = Cj for all j G a-i and z,; • a,- < Cj for 
all j $l o~i. The sets a\, . . . , at are called the cells of the regular subdivision and the regular- 
subdivision is denoted by A C (A) = {a±, . . . , at} or simply A c when A is unambiguous. 

Equivalently, regular subdivisions are simply capturing complementary slackness from 
linear programming. Namely, a feasible solution to LP,4 jC (b) is optimal if and only if the 
support of the feasible solution is a subset of some cell of A c . Geometrically, A c can be 
thought of as a partition of cone(^4) by the inclusionwise maximal ones among the cones 
cone ( A ai ),..., cone (A at ); each such cone is generated by the normal vectors of defining 
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inequalities of faces of Q c , each maximal cell indexes the set of normal vectors of inequalities 
satisfied with equality by a vertex (or minimal face) of Q c . The zi, . . . ,Zf above are these 
vertices and so the regular subdivision A c is geometrically realized as the normal fan of Q c . 

A regular subdivision of A is called a triangulation if the columns of each A a . are linearly 
independent for alH = 1, . . . , t. Note that a regular subdivision A c is a triangulation if and 
only if every vertex is contained in exactly d facets; that is, the polyhedron Q c is simple, 
or, non-degenerate. A triangulation A c is called unimodular if det{A ai ) = ±1 for each 
maximal cell of A c . A refinement of a subdivision A c of A is another subdivision A c / of 
A so that each cell of A c / is contained in some cell of A c . A set B C Z d is a Hilbert basis 
if N-B = cone(-B) n 1j d ; we will say that a matrix B is a Hilbert basis if its columns form 
a Hilbert basis. Note that if for some c £ W 1 A c is a unimodular triangulation of A then 
Cramer's rule implies that A itself is a Hilbert basis. 

Let IPa,c := {IPA,c(b) : b £ ~NA} denote the family of integer programs IP^^b) 
having a feasible solution. Informally, a test-set for the family of integer programs IPac * s a 
finite collection of integer vectors, called test vectors, with the property that any non-optimal 
feasible solution to any integer program in this class can be improved (in objective value) 
by subtracting a test vector from it. Test-sets for the family IPa,c were first introduced 
by Graver [18]. Graver's test set, called a Graver basis in the literature, is typically not a 
minimal test set. 

A simple but helpful characterization of the TDI property in terms of the Hilbert basis 
property of regular subdivisions has been provided by Schrijver [33]. We prove another 
elementary characterization in Section [2] whose simplified version is the following: 

Theorem 12.51 The system yA < c is TDI if and only if A is a Hilbert basis, and there 
exists a test-set for IPa,c with all test vectors having positive entries equal to 1 and this 
positive support indexes either a linearly independent set or a minimally linear dependent 
set of columns of A. 

Establishing this theorem, some of its corollaries (see Section [2]) have led us to three 
relevant earlier results that have been found in contexts different from ours, independently 
of one another: 

- Applegate, Cook &; McCormick's result [H Theorem 2] which states that if A is a 
Hilbert basis then yA < c is TDI if LP,4 jC (b) has an integer optimum for every 
constraint vector b that is a sum of linearly independent columns of A. However, 
our test-sets and their use for optimizing on IP are new. 

- Hosten and Sturmfels result |21[ Theorem 1.1] is stated in a computer algebra con- 
text, but the special case when the gap is can be seen to be equivalent to the 
above statement of Theorem 12.51 However, it is still easier to prove each of these 
two independent results in their own terms, directly. 

- Eisenbrand and Shmonin [13] proved Ho§ten and Sturmfels result in an elementary 
way, yet in the context of a generalization of integer programming called "parametric 
integer programming". However, this work also uses more involved machinery than 
we do here. 
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In addition, none of these works characterizes the TDI property in terms of test-sets or 
solving the underlying integer programs in fixed dimension without using Lenstra's algo- 
rithm. As previously noted in [T] , this allows us to deduce shortly Cook, Lovasz & Schrijver's 
result [7J on testing for the TDI property in fixed dimension, when given the conditions of 
Theorem [231 

Another virtue of Theorem 12.51 is a useful reformulation to polynomial ideals: it general- 
izes a well-known algebraic result proved by Sturmfels [381 Corollary 8.9] relating toric initial 
ideals to unimodular triangulations. The basic connections between integer programming 
and computational algebra, knowledge of which will not be assumed here, was initiated by 
Conti and Traverso [5] and further studied from various viewpoints in [37], [38], [39], [40] . 

If A is a matrix whose first d x (n — d) submatrix is a — 1 matrix and whose last d x d 
submatrix is —Id, and c is all 1 except for the last d coordinates which are 0, then DP/i )C (b) 
is called a set packing problem, and Q c a set packing polytope. In Section [3] we show that the 
converse of the following fact (explained at the end of Section [2j holds for normal fans of 
integral set packing polytopes: if c,c' 6 MJ 1 are such that A c i is a refinement of A C; where 
A c / is a unimodular triangulation, then yA < c is TDI. In general, the converse does not 
hold and the most that is known in this direction is the existence of just one full dimensional 
subset of the columns of A which is unimodular [16]. Not even a "unimodular covering" 
of a Hilbert basis may be possible [2j. However, the converse does hold for normal fans of 
integral set packing polytopes. More precisely, the main result of Section [3] is the following: 

Theorem 13.11 Given a set-packing problem defined by A and c, Q c has integer vertices if 
and only if there exists c' such that A c / is a refinement of the normal fan A c of Q c , where 
A c / is a unimodular triangulation. 

The proof relies on the basic idea of Fulkerson's famous "pluperfect graph theorem" [17] 
stating that the integrality of such polyhedra implies their total dual integrality in a very 
simple "greedy" way. Chandrasekaran and Tamir [3] and Cook, Fonlupt and Schrijver [6] ex- 
ploited Fulkerson's method by pointing out its lexicographic or advantageous Caratheodory 
feature. In [35\ §4] it is noticed with the same method that the active rows of the dual of 
integral set packing polyhedra (the cells of their normal fan) have a unimodular subdivision, 
which can be rephrased as follows: the normal fan of integral set packing polyhedra has a 
unimodular refinement. However, the proof of the regularity of such a refinement appears 
for the first time in the present work. 

These results offer four methods for recognizing TDI systems, explained, illustrated and 
compared to previously known ones as well, in Section HI Particular attention will be given 
to the recognition of integral set packing polytopes in practice, through Theorem 13. II above, 
and by using computational algebra packages like Macaulay 2. 

2. TDI Systems 

In this section we provide some new characterizations of TDI systems. We show the 
equivalence of five properties, three polyhedral (one of them is the TDI property) and two 
concern polynomial ideals. A third property is also equivalent to these in the generic case. 
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While the proofs of the equivalences of the three polyhedral properties use merely poly- 
hedral arguments, the last among them - (iii) - has an appealing reformulation into the 
language of polynomial ideals. Therefore, we start this section by introducing the necessary 
background on polynomial ideals; namely, toric ideals, their initial ideals and Grdbner bases. 
The characterizations of TDI systems involving polynomial ideals are useful generalizations 
of known results in computational algebra. See [TT] and [38] for further background. 

An ideal I in a polynomial ring R := \c[xi, . . . ,x n ] is an R- vector subspace with the 
property that I ■ R = I. It was proven by Hilbert that every ideal is finitely generated. 
That is, given an ideal I there exists a finite set of polynomials fx, . . . , ft G I such that for 
every / G I there exists hi, . . . , hf G R with / = h\f\ + • ■ • + h t ft- We call such a collection 
fi, . . . , ft G la generating set for the ideal / and denote this by I = (fi, . . . , ft). For the 
monomials in R we write x u = x" 1 • • • x™" for the sake of brevity. We call u the exponent 
vector of x u . A monomial x u is said to be square-free if u G {0, l} n . An ideal is called a 
monomial ideal if it has a generating set consisting only of monomials. For any ideal J of 
R, mono( J) denotes the monomial ideal in R generated by the set of monomials in J. An 
algorithm for computing the generators of the monomial ideal mono (J) can be found in [321 
Algorithm 4.4.2]. 

Every weight vector c G ffi" induces a partial order y on the monomials in R via x u y x v 
if c • u > c • v. If c G M n where 1 is the monomial of minimum c-cost (that is, c • u > for 
every monomial x u ), then we can define initial terms and initial ideals. Given a polynomial 
/ = X^ugN" r uX u G I the initial term of / with respect to c, is denoted by in c (/), and 
equals the sum of all r u x u of /, where c • u is maximum. Note that we can always write 
a polynomial / as / = in c (/) + trail c (/) in the obvious way. The initial ideal of I with 
respect to c is defined as the ideal in R generated by the initial terms of the polynomials 
in /: in c (/) := (in c (/) : / G /). A Grdbner basis of an ideal I with respect to c, is a 
finite collection of elements g±, . . . ,g s in I such that in c (I) = ( in c (gi), in c (g2), • • • , in c (ff s ) )• 
Every Grobner basis is a generating set for the ideal I. 

If in c (I) is a monomial ideal then a Grobner basis is reduced if for every i ^ j, no term 
of gi is divisible by m c (gj). The reduced Grobner basis is unique. In this case, the set of 
monomials in in c (I) equal {x u : u G U} with U := D + N™ where D is the set of exponent 
vectors of the monomials in c (g , i), in c (g2), • • • ,^c(Ss)- Dickson's lemma states that sets of 
the form Z) + N n , where D is arbitrary have only a finite number of minimal elements (with 
respect to coordinate wise inequalities). This is an alternative proof to Hilbert 's result that 
every polynomial ideal is finitely generated. In this case, the Grobner basis also provides a 
generalization of the Euclidean algorithm for polynomial rings with two or more variables 
called Buchberger's algorithm (see Chapter 2, §7]). This algorithm solves the ideal 
membership problem: decide if a given polynomial is in an ideal or not. However, a Grobner 
basis for an ideal can have many elements (compared to a minimal generating set for the 
ideal), and none of the related computations can be achieved in polynomial time. 

If in c (I) is not a monomial ideal then we can form a monomial initial ideal, in^ c (I), as 
follows: fix an arbitrary term order independent of c, (that is, a total ordering >- of the 
vectors in N" - satisfying u y for every u G N" - , and if u y v then u + 7 y v + 7 for all 
7 G N n ). We use this term order to break ties: x u ^ c x v if and only if c • u > c • v, or 
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c • u = c • v and u >- v. Clearly, 'V c " is also a term order and so in>_ c (/) is a monomial 
ideal. Such a tie-breaking will be needed in the proof of Lemma 12.21 

The toric ideal of A is the ideal I a = (x u — x v : Au = Av, u, v € N n ) and is called 
a binomial ideal since it is generated by polynomials having at most terms. Every reduced 
Grobner basis of a toric ideal consists of binomials. A toric initial ideal is any initial ideal 
of a toric ideal. 

Remark 2.1. It follows for (a not necessarily generic) c, that the reduced Grobner basis of 
a toric ideal I a (with respect to >- c ) is of the form Qy c = {x u ? — x u ^ : i = 1, . . . , t} and we 
can suppose that in c (x u ^ — x u * ) = x u ^ for i = 1, . . . , s and in c (x u ^ — x u * ) = x u ^ — x u * 
for % = s + 1, . . . ,t. Furthermore, the set of polynomials Sy := {x u ^ , . . . ,x u ^,x Us + 1 — 
x Us +! , . . . , x u ^ — x u * } is a reduced Grobner basis for in c (/^) with respect to the term order 
>- cf. [68\ Corollary 1.9]. Note that if c were generic then s = t and Sy is simply the 
minimal generating set for in^i^). 

The following lemma is a natural connection between integer programming and toric 
initial ideals. It originally appeared in [32\ Lemma 4.4.7] but we prove it here in order for 
this article to be self-contained. 

Lemma 2.2. For A G Z dxn and c G M. n the monomial ideal mono(in c (/ J 4)) is equal to 
(x w : uj G N n is a non-optimal solution for IP a, c (Alu) ). 

Proof. It is straightforward to show that mono(in c (I^)) contains the defined set: let u be 
a non-optimal solution, and u' an optimal solution to IPa,c(Aui). Then x w — x w G I a is a 
binomial having the monomial x w G mono(in c (Iyi)) as its initial term with respect to c. 

Suppose now x^ G in c (/^). That is, there exists a polynomial / G Ia with in c (/) = x w . 
We will prove that u is a non-optimal solution for IP^^A^). 

Let >- be an arbitrary term order and let Sy be the reduced Grobner basis of in c (I^) 
with respect to >- as in Remark 12. 1[ We proceed by induction with respect to the minimum 
number of successive polynomial divisions of / by the elements of Sy, and replacing / by 
the remainder after each division, until arriving at a remainder. As noted in Remark 12.11 
Sy = {x< , . . . , x u ^ , x u ^+i - x u s~+i , . . . , x u ^ - x u *~ }. 

If first we divide / with a monomial in Sy , say x u i^ , then c • > c • , and in c (/) = x w 
is divisible by x u ^ , so uj — uf > 0, c- (u — (uf — uj~)) < c-uj, where A{ui — (uf — u^)) = Alo. 
Threrefore uj is also a non-optimal solution to IPa,c(Auj). 

Otherwise, the first polynomial division is by a binomial in Sy, say, x u ^ — x Ut . Now 
recall that / G I a is a polynomial with in c (/) = in^ c (/) and so the resulting polynomial 
after the division can be written as /' := -^rx 11 ' — trail c (/). The initial term in c (/') is a 

monomial and equals ^^x u « since c • (u — (uf — u^T)) equals c ■ uj, due to c • uf = c • u^. 

This implies that in c (trail c (/)) is strictly cheaper (with respect to c) than ^rx u * and 
so -^x"' G in c (I^). Since /' requires one less division than / by Sy to arrive at 

-, u + 
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remainder then, by induction, we have that uj — (uf — u t ) is not an optimal solution to 
IPa,c(^4(w — ( u ^ — u <~))) an d so w i s n0 ^ an optimal solution to IPa,c(Au>) either. □ 

Remark 2.3. a C [n] is contained in a cell of A c if and only if ^ ig(T e.j E R n is an optimal 
solution to LPyi iC (b -) where b CT := ^4(X^ gCT e i)- This happens in particular if a = 0. Indeed, 
then Ei GCT e i = °) and b <r := -4(Ei G(T e i) = °- For this b = G K d , £ R" is an optimal 
solution, otherwise LP J 4 c (b) is not bounded for any b, and Q c = 0. We are not interested 
in such polyhedra, and we will avoid this situation by supposing Q c / 0. 

A test-set [18] for the family of integer programs IP^.c := {IPA,c(b) : b E NA} is a 
collection of integer vectors {v^~ — : Avf = Av~ , c • v; + > c • Vi~, vf, v 4 ~ E N n , i = 
1, . . . , s} with the property that for all feasible, non-optimal solution u to IP^ ;C (b) there 
exists an i, 1 < i < s, such that u — (v^~ — v~) > 0. We can now state and prove our 
characterizations of TDI. 

Example 2.4. Let us suppose Q c ^ 0, that is, E M n is an optimal solution for LP^ jC (0), 
and show a particular test-set related to Remark 12.31 We will say that k C [n] is a wheel 
(with respect to A and c), if A(J2 i&K ei) = 0, X^«e«; Ci > ^ ano - ^ or all i € k, k \ {i} is a 
subset of a cell (of A c (^4)). 

For instance if A = (ai, . . . , a^, — (ai+. . .-t-a^)), where linearly independent 

integer vectors, and c > 0, then A c (^4) is a triangulation whose maximal cells are precisely 
{[d + l]\{i} : i = 1, 2 . . . d + 1}. In this case the one-element set {(1, . . . , 1) E R d+l } is a 
test-set for IPa,c- Note that if cone(^4) is pointed then A has no wheel, regardless of the c 
in question. 

Theorem 2.5. Fix A E Z a!xn , where A is a Hilbert basis, c E M n , and Q c / 0. The 
following statements are equivalent: 

(i) The system yA < c is TDI. 

(ii) T/ie subconfiguration A a of A is a Hilbert basis for every cell a in A c . 

(iii) There exists a test- set for TP a,c where all the positive coordinates are equal to 1, and 
each positive support is either the incidence vector of a linearly independent set of 
columns, or of a wheel; in the former case the negative support is a subset of a cell, 
in the latter case it is 0. 

(iv) The monomial ideal (x^ : u £ N n is not an optimal solution for IPa,c(Auj) ) has a 
square-free generating set. 

(v) The monomial ideal generated by the set of monomials in in c (/ J 4), that is, mono(in c (/yi)) 
has a square-free generating set. 

The main content of the equivalence of (ii) and (v) is that the "Hilbert basis property" 
is equivalent to the existence of a "square-free generating set", extending a well-known 
result |38|. Corollary 8.9] (see Corollary 12.71 below) to the non-generic case. We have re- 
cently found a similar extension in Hosten and Sturmfels [21], equivalent to ours. (The 
latter paper provides an algorithm to compute the integer programming gap gap^ c := 
max{OPTLP Ac (b) - OPTLP Ac (b) : b £ NA}, where OPTIP and OPTLP mean the op- 
timal value of the corresponding programs. Note that the system yA < c being TDI is 
equivalent to A being a Hilbert basis and gap^ c = 0. Ho§ten and Sturmfels compute 
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§ a PA c by studying the "primary decomposition" of the monomial ideal mono(in c (/ y 4)). We 
present below our direct elementary argument.) 

Note also that (i) implies that A is a Hilbert basis, since this latter is equivalent to the 
following property: if LP^ c (b) has a feasible solution, then it also has an integer feasible 
solution; (ii) also implies it, since a cone that has a subdivision to Hilbert cones is itself a 
Hilbert basis. So this condition has to be added only to (iii), (iv) and (v); they may all be 
satisfied without A being a Hilbert basis, certifying then that yA < c is not TDI, unless 
Q c = 0. The condition Q c ^ makes sure that LP J 4. c (b) is bounded, otherwise already the 
statement of (iii), (iv), (v) is not well defined. 

Proof, (i) => (ii) : This is well-known from Schrijver's work, (see for instance [33])) but 
we provide the (very simple) proof here for the sake of completeness: Suppose the system 
yA < c is TDI, and let a € A c . We show that A a is a Hilbert basis. Let b 6 cone(A a ). 
Since the optimal solutions for LP^ jC (b) are exactly the non- negative combinations of the 
columns of A a with result b, the TDI property means exactly that b can also be written 
as a non-negative integer combination of columns in A a , as claimed. 

(ii) =4* (iii) : Suppose (ii) holds true for A c of A. For every r C [n] with r not contained in 
any cell of A c , let b T := Y2ie T a i = A(J2ier e i)- Since r is not contained in any cell of A c , 
there exists an optimal solution (5 T to LP^ >c (b r ) with c • j3 T < c • Ylier e «- By the optimality 
of (3 T we must have supp(/3 r ) C a for some cell a of A c (^4). 

If b r ^ 0, then a ^ 0, and by (ii) A a is a Hilbert basis. Therefore (3 T can be chosen to 
be an integral vector. On the other hand, if b r = 0, then by the condition Q c ^ 0, is an 
optimal solution for LP J 4 jC (b r ) (see Remark 12. 3ft . Let 

Ta,c := { e,j — [3 T : A T is linearly independent and r is not contained in any cell of A c } 
[^J ■ k is a wheel (with respect to A and c)}. 

We claim that Ta,c is a test-set for IPy^c- By construction every t € Ta,c satisfies At = 
and c • t > 0, so we have to prove only the following: 

Claim: For every b G Z d and feasible but not optimal solution lo o/IPyi,c(b), there exists 
t G T~A,c such that lo — t > 0. 

We have Aoj = b, where supp(w) is not contained in any cell a of A c . 
Case 1. b := Aw ± 0: 

By basic linear programming ("Caratheodory's theorem") there exists r C supp(oj), 
so that A T is linearly independent, and b £ cone(^4 T ). If every r satisfying these three 
conditions is a subset of a cell, then the unique solution of the equation A T x T = b is feasible 
for LPyi )C (b) and optimal, and again, by basic linear programming this would contradict 
that u is not optimal. So r can be chosen so that it satisfies these three properties, but it is 
not contained in any cell o/A c . Clearly, uj — (X^g T e i~ftr) > 0, and t := X^er e i~ftr G Ta,c-, 
finishing the proof of the claim. 

Case 2. b := Au = 0: 
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Since then is optimal, and w is not, we have 

c-w>0. (1) 

(a) every proper subset o/supp(cj) is contained in some cell, and w is a — 1 vector. 

First, if there exists w' G N n , w' <w, with u/ not optimal for / b' := Aw', then we can 
apply Case 1 to w', and get t G Tj± c , u — t > u>' — t > 0, and we are done. Second, if w has 
a coordinate that is bigger than 1, say w± > 1, then Case 1 can be applied to a/ := a; — ei . 

We can now assume that (a) holds. Letting k := supp(u;), we show: 

(b) every proper subset of A K is linearly independent. 

Suppose (b) is not true. Then for some J C k there exists a coefficient vector A = 
(Ai, . . . , A n ) G M n , Xj ^ if j G J, and Aj = if j £ J. Moreover (by dividing with the 
highest coefficient, say that of the first coordinate), Aj < 1 (J G [n]), and Ai = 1, such that 

AX = 0. 

Furthermore, we can partition J = J + U J~ where J + := {j : Aj > 0} and J~ := {j : 
Aj < 0}. Using this partition, we can rewrite A in the obvious way as A = A + — A~ where 
A + , A - > with the supports being J + and J~ respectively. Consequently, we have 

^A + = .4A~ 

and according to (a) there exist cells of A C (A) containing J + and J~ so both A + and A~ 
are optimal solutions to the same linear program implying that: 

c.A + = c-A~. (2) 

By our choice of the size of the components in A = A + — A™, we have Ai = 1. In addition, 
A < w and w\ = 1 and, since supp(o; — (A + — A~)) C C k then, by (a) and the fact 

that AX + = AX~, w - (A + - A") is an optimal solution to LPy!i )C (b). 

Recall that is also an optimal solution to LPyi iC (b). Combining (1) and (2) above, we 
get that c ■ (w — (A + — A - )) > = c • which contradicts the optimality of w — (A + — A - ) 
and so no such linear dependent JCr can exist. This proves (b). 

Now by (a) and (b): k is a wheel for A and c, and as such ^igK e « ^ Ta,o finishing the 
proof of the claim and of (ii) => (iii). 

(iii) =>■ (i) : Suppose (iii) holds, and b G cone(A). Then LP J 4 iC (b) is feasible, and since 
Q c ^ the optimum is bounded. So by the condition that A is a Hilbert basis, IP^ ^b) 
is also feasible and it also has a bounded minimum attained by w G N n . Suppose for a 
contradiction that the optimal solution a/D to LPyi iC (b), (a G N™, D is a positive integer) 
satisfies c • a/D < c • u. This also implies that Dw is not an optimal solution to IP^ c (J5b). 

By (iii) there exists a 7 + — 7~ G Ta,c with 7 + G {0, 1}™ and 7" G N n such that 
c • (7+ — 7~) > and Dw — (7+ — 7") G N n . Hence, supp(7 + ) C supp(Z)u;) = supp(u;). 
Since the value of all elements in 7 + is or 1 then we also have w > 7 + , so w— (7 + — 7~) G W 1 
is also a feasible solution to IPA, c (b) with c • (w — (7 + — 7")) < c • w, in contradiction to 
the optimality of w. 
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(iii) 44> (iv) 44> (v): Both (iii) and (iv) can be reformulated as follows: If uj 6 N n is not 
an optimal solution to IVa,c{Auj) then the vector u>' := X^iesupp(w) e * * s a ^ so a non-optimal 
solution to IPa,c(Auj'). The equivalence of (iv) and (v) is a special case of Lemma l2.2i □ 

The simple equivalence of (i) with (iii) is a possible alternative to Cook, Lovasz &; 
Schrijver's result [7J, and (iii) could also be replaced by a characterization of Applegate, 
Cook & McCormick pQ stating that it is sufficient to check the property for functions b 
that are small subsums of rows of A. For the sake of completeness and for later reference 
we present the corresponding computational corollary in our context: 

Corollary 2.6. [T], [7J Let the dimension d be fixed but the system yA < c, given as input, 
have an arbitrary number of inequalities, where A is a Hilbert basis. Then this system can 
be tested for the TDI property in polynomial time. 

Proof. We use the equivalence of (i) and (iii) and the construction of Ta.c from above. Note 
that the wheels are straightforward to identify but they do not have to be identified for 
testing TDI-ness. 

We can construct the non-wheels of Ta,c or conclude that the system is not TDI in 
0(n d ) time. Listing all the linearly independent (or generously all the d element) subsets r 
of {1, . . . , n} that are not subsets of cells in the fan. For each such subset we can form the 
vector b T := ^4(X^ei- e «) ana - identify its cell a in the fan by solving LPyi iC (b r ). 

Next, for each such r, with its corresponding optimal cell a, do the following: we can 
either find an integer optimal solution (3 T to LPyi jC (b T ) (which will have support in a), or 
else no such integer solution exists in which case we conclude that yA < c is not TDI, which 
can be done in polynomial time - see [20|. §6.7] for example. 

If LP^ jC (b T ) has an integer solution then repeat the same test, as for b T above, for 
b r — rjaj where i £ a and is the largest positive integer for which the optimal solution to 
LPyi jC (b r — rjaj) is in the same cell as a. If LP J 4 iC (b r — r^i) has an integer optimum solution 
then repeat as above by replacing b r with b r — rjaj. If at any point during this procedure 
the updated b T yields no integer optimal solution to LP J 4 jC (b r ) then we immediately stop 
and conclude that yA < c is not TDI. Otherwise, we add ^2 iGr ^% — T to the test set. 

We repeat this procedure for every such r as above and if for none of the linearly 
independent sets of columns r (different from the cells) we arrive at the conclusion that the 
system is not TDI, then we have in 7^ jC an element of the form ^2 i€r ej — (3 T for every r 
that is not subset of a cell, and we can conclude like in the proof of the theorem, that Ta,c 
is a test-set. Therefore (iii) is satisfied, so according to Theorem 12.51 (i) is also satisfied, 
that is, yA < c is TDI. If d is fixed, A has a polynomial number of linearly independent 
subsets of columns. □ 

One consequence of this corollary is that the constructed test-set allows us to solve 
IP.4,c(b) in polynomial time without the use of Lenstra's algorithm. 

Recall that we defined c £ W 1 to be generic with the first of the following conditions; 
the others are equivalent properties for toric ideals [40j : 

- The integer program IP^ ;C (b) has a unique optimal solution for all b S NA. 
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- The toric initial ideal in c (lA) is a monomial ideal. 

- There exists a reduced Grdbner basis {x Ul — x u i , . . . , x Us — x Us } of I a with c • > 
c • for each i = 1, . . . , s. 

Recall that if each IP^.^b) has a unique solution then so does LP^ jC (b). Hence in 
the generic case, assuming Q c ^ 0, Cramer's rule tells us that A c being a unimodular 
triangulation of A implies that A is a Hilbert basis without supposing it in advance; and so 
(ii) implies (v) without any condition; to prove the converse statement in the generic case, 
(v) along with the assumption that 'LA = 7j d implies that A is a Hilbert basis, and then we 
can apply again Theorem 12.51 (v) implies (ii) to get: 

Corollary 2.7. (Sturmfels) [38j Corollary 8.9] Let A £ Z dxn and let c £ R n be generic with 
respect to A. Then A c is a unimodular triangulation if and only if the toric initial ideal 
in c (/^) is generated by square-free monomials. 

Theorem [23] is the result of generalizing Sturmfels' above result to arbitrary TDI systems. 
Still concerning generic c it is worth to note the following result of Conti and Traverso which 
provides another connection between integer linear programming and Grobner bases. Here 
we think of an element x v — x v as a vector v + — v — . 

Proposition 2.8. (Conti- Traverso) [3] - see [HJ Lemma 3] // IP,4,c(b) has a unique 
optimal solution for every b £ N^4 then the reduced Grobner basis is a minimal test-set for 
the family of integer programs IPa,c- 

This proposition means for us that in the generic case the following (vi) can be added 
to Theorem 12.51 (vi) The initial terms in the reduced Grobner basis are square-free. In 
particular, in the generic case of condition (hi) of Theorem 12.51 the unique inclusionwise 
minimal test-set is defined by the reduced Grobner basis, which, by (vi) has only square- 
free terms initial terms. Even though Theorem 12.51 concerns general TDI systems and could 
be proved in elementary means, it was highly stimulated by the above results concerning 
the generic case. 

As is typically the case in combinatorial optimization, the cost vector c is not generic 
for A. However, there may be cases where one can slightly perturb c to another cost vector 
that is generic for A, and where the TDI property for yA < c can be more easily studied 
when c is perturbed. More precisely, from the implication "(ii) implies (i)" we immediately 
get the following: 

Proposition 2.9. // c, c' £ IR n are such that A c / of A is a refinement of A c of A, where 
A a is a Hilbert basis for all a £ A c /, and in particular if A c / is a unimodular triangulation 
of A, then yA < c is TDI. 

Remark 2.10. Having a regular unimodular refinement A c / of A c amounts to providing 
an integer point on the face of each P^, for all b £ cone(^4) n that is minimized by c. 
This integer point is the vertex of Pb minimized by c' and so having a regular unimodular 
refinement provides an integer point on each Pb in a uniform manner, dictated by c' . 

In the rest of the paper one of our favorite themes will be the use of Proposition 12.91 
Clearly, the unimodular triangulation does not even need to be regular: in fact, a unimodular 
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cover of each of the cells of A c suffices for verifying that yA < c is TDI. However, we are 
interested in cases when the converse of Proposition 12,91 is true. In general it is not. It is 
not even true that a Hilbert basis has a unimodular partition or a unimodular covering [2] 
and this counterexample inspires two more remarks. First, it cannot be expected that the 
equivalence of (i) and (v) can be reduced to Sturmfels' generic case. Secondly, it should be 
appreciated that the converse of this remark does hold in the important set packing special 
case, as we will see in the next section. 

3. Set Packing 

Let a set packing problem be defined with a matrix A and vector c, and recall c := 
(1,0) € R n , where the last d entries of c are 0. If the set packing polytope Q c has integer 
vertices then the matrix A and the polytope Q c are said to be perfect. (We will not use the 
well-known equivalence of this definition with the integer values of optima: this will follow.) 
Lovasz' (weak) perfect graph theorem [25] is equivalent to: the matrix A defining a set 
packing polytope is perfect if and only if its first (n — d) columns form the incidence vectors 
(indexed by the vertices) of the inclusionwise maximal complete subgraphs of a perfect graph. 

A polyhedral proof of the perfect graph theorem can be split into two parts: Lovasz' 
replication lemma [25] and Fulkerson's pluperfect graph theorem |17j . The latter states 
roughly that a set packing polytope with integer vertices is described by a TDI system of 
linear inequalities. In this section we restate Fulkerson's result in a sharper form: there is 
a unimodular regular triangulation that refines the normal fan of any integral set packing 
polytope. We essentially repeat Fulkerson's proof, completing it with a part that shows 
unimodularity along the lines of the proof of [35L Theorem 3.1]. The following theorem 
contains the weak perfect graph theorem and endows it with an additional geometric feature. 

Theorem 3.1. Let Q c be a set packing polytope defined by A and c. Then there exists a 
vector e € W 1 such that c' := (1, 0) + e defines a regular triangulation A c / refining A c , and 
this triangulation is unimodular if and only if Q c is perfect. 

We do not claim that the following proof of this theorem is novel. All essential ingredi- 
ents except unimodularity are already included in the proof of Fulkerson's pluperfect graph 
theorem [17] . Fulkerson's proof suggests a greedy way of taking active rows with an integer 
coefficient (see below); this is often exploited to prove that some particular systems are 
TDI, let us only cite two papers the closest to ours, Cook, Fonlupt & Schrijver [6] and 
Chandrasekaran & Tamir [3j. The latter paper extensively used lexicographically best so- 
lutions, which is an important tool in linear programming theory. This idea was used in 
|35j to prove that the normal fan has a unimodular refinement. This same lexicographic 
perturbation is accounted for by the vector e of Theorem 13. 1[ showing that the unimodular 
refinement is regular. This motivated the following problem, which contains the perfectness 
test (detecting perfection in a graph): 

Problem 3.2. |36j Given a d x n integer matrix A and an n dimensional integer vector 
c, decide in polynomial time whether the normal fan of Q c consists only of unimodular 
cones. Equivalently, can it be decided in polynomial time that Q c is non- degenerate, and 
the determinant of A a is ±1 for all a E A c . 
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Motivated by the perfectness test, the following problem might still be polynomially 
solvable. 

Problem 3.3. Given a d X n integer matrix A and an n dimensional integer vector c, 
decide in polynomial time whether A c has a unimodular refinement A c / . 

In the set packing case not only is the perfectness test more efficient, but also the linear 
progams and their duals can be solved in polynomial time (even if the algorithm uses 
nonlinear optimization and the ellipsoid method). Could this also be true in general ? 

Problem 3.4. Given a d x n integer matrix A and an n dimensional integer vector c as 
input, is there a polynomial algorithm that finds at least one of the following as output: a 
solution to lP^ c (h), or a "NO" answer to Problem \3.3\ (or to Problem \3.2\) . 

All these problems can be solved in polynomial time in the set packing case: Problem l3.3l 
according to Theorem [3J] for the moment only by using [9]; Problem l3.4l for the moment only 
by [2D] . It remains an interesting question whether there are more efficient and conceptually 
simpler solutions to these problems. 

We now prepare the proof of Theorem 13.11 It is a last step in a sharpening series of 
observations all having essentially the same proof. We begin with the proof of Fulkerson's 
pluperfect graph theorem which will indicate what the c' of Theorem 13.11 should be, and 
then finish by showing that A c / is a unimodular triangulation. 

Assume that A is a perfect matrix for the remainder of this section and that c = (1, 0) 
as before. For all b £ Z rf and i £ [n] let 

A c i(b) := maxjxj : x is an optimal solution of LP^ c (b) }. 

That is, A Cj j(b) is the largest value of Xi such that c • x is minimum under x £ P^. 

Remark 3.5. If a is the minimal cell of A c such a b £ cone( J 4 (J ), then b — A C] j(b)aj £ 
cone( J 4 (J ') where a 1 £ A c , a' C a and the dimension of cone(A a /) is strictly smaller than 
that of cone(A a ). Furthermore, b — Aaj £" cone(A a ) if A > A c/ ;(b). 

For all b £ Z d we show that A Cj j(b) is an integer for every i = 1, . . . , n. This is the heart 
of Fulkerson's pluperfect graph theorem [17} Theorem 4.1]. We state this in the following 
lemma in a way that is most useful for our needs. Denote the common optimal value of 
hpA,c(b) and DP J 4 ;C (b) by 7 c (b). Note that j c is a monotone increasing function in all of 
the coordinates. 

Lemma 3.6. Suppose 7 c (b) £ Z for all b £ Z d . // x is an optimal solution to LP J 4 )C (b) 
with X[ ^ for some 1 < I < n, then there exists x* also optimal for LP J 4 iC (b) ; such that 
x* > 1. 

Note that this lemma implies the integrality of A := A Cj /(b) for all I = 1, . . . , n: if A 
were not an integer then setting b' := b — [Aja/ we have A Cj /(b') = {A} where < {A} := 
A — [AJ < 1 , contradicting Lemma 13.61 

Proof. Suppose x £ P^ with c x = 7(b) and xi > for some 1 < I < n. We can assume that 
X[ < 1 since otherwise x* := x. We have two cases: either 1 < I < n — d or n — d+1 < I < n. 
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If n — d + 1 < I < n then = — e;_( n _^ G M d and q = 0. In this case, we have 
7c(b) = 7c(b + xiei_(, n _ d ^) because replacing x\ by in x we get a solution of the same 
objective value for the right hand side b + x;e;_( n _ rf ) which gives 7 c (b) > 7 c (b + x/e/„( n _ rf )). 
The reverse inequality follows from the (coordinate- wise) monotonicity of j c . But then 

7c(b + e Mn _ d) ) < 7 c (b + :qe Mn _ d) ) + 1 - x t = 7 c (b) + 1 - x h 

and since 7 c (b + e^_( n _ rf )) is integer and < 1 — X\ < 1, we conclude that 7c(b + e i „( n _ d )) = 
7c(b). 

So for any optimal x' G Pb+e,_ (n _ d) where c • x' = 7 c (b), letting x* := x' + ei-(n~d) G Pb 
we have c • x* < 7 c (b) and so x* is optimal and x^ > 1. 

Suppose now 1 < I < n — d. Replacing xi in x by we get a point in Ph-x^ ■ This point 
has objective value c • x — xi < c • x = 7 c (b), and so we have by monotonicity 

7(b - a z ) < 7(b - xiaf) < 7(b). 

Since the left and right hand sides are both integer values then 7(b — a;) < 7(b) — 1. Letting 
x* := x' + ei G Pb we have c • x* < 7 c (b) — 1 + 1 = 7 c (b), so x* is optimal, and x^ > 1. □ 

Let us now define the appropriate c' for the theorem, depending only on c. Define 
c' := c + e G W 1 where £j := — {l/d d+2 ) 1 for each i = 1, . . . , n. Note that the absolute value 
of the determinant of a d x d {—1,0, l}-matrix cannot exceed d d . It follows, by Cramer's 
rule, that the coefficients of linear dependencies between the columns of A are at most d d in 
absolute value, and then the sum of absolute values of the coefficients between two solutions 
of an equation ^4x = b for any b G M. d can differ by at most a factor of d d+2 . After this 
observation, the following lemma is straightforward to verify 

Lemma 3.7. (i) Any optimal solution to LP J 4 c /(b) is also optimal for LP y4,c(b) . 

(ii) // x' and x" are both optimal solutions to LP,4 )C (b) then x' is lexicographically 
bigger than x" ( that is, the first non-zero coordinate of x' — x" is positive ) if and 
only if c' • x' < c' • x" . 

Statement (i) of the lemma means that A c / refines A c , and statement (ii) means that 
an optimal solution to LP^ c ^(b) is constructed by defining b° := b and recursively 

Xi := A Ci j(b i_1 ), b* := b* _1 — for alii = 1, . . . ,n. 

Furthermore, this optimum is unique and it follows that A c / is a triangulation. We are now 
ready to prove Theorem 13. 11 

Proof of Theorem \'J.l\ The necessity of the condition is straightforward: each vertex y G Q c 
satisfies the linear system (consisting of d equations in d unknowns) of the form yA a i = c a i 
where a' is a cell of A c /, c CT / is the subvector of c indexed by a'. Since the determinant of 
A a r is ±1, y must be an integer vector because of Cramer's rule. 

Conversely, we will prove the assertion supposing only that 7 c (b) is integer for all b G Z d 
and applying Lemma 13.61 Note that, by the already proven easy direction, we will have 
proved from this weaker statement that Q c is perfect. Without loss of generality, suppose 
that b G 1j d cannot be generated by less than d columns of A. That is, the minimal cell a 
of A c such that b G cone(^4 CT ) is a maximal cell of A c . That is, cone{A a ) is d-dimensional. 
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Because of Lemma l3.7f i). an optimal solution to LPyi jC '(b) will have support in a and 
Lemma l3.7f ii) implies that such an optimal solution is constructed as follows: Let s\ := 
min{i : i G a} and x si := A C)S1 (b). Recursively, for j = 2,...,d let Sj be the smallest 
element in a indexing a column of A on the minimal face of cone(j4<j) containing 

i-i 
i=l 

Since b is in the interior of cone(A CT ) then x Si > for each i = l,...,d, and by 
Lemma |3.6[ these d x Si 's are integer. Moreover, since the dimension of cone(A CT \ r si ai \ ) 
is strictly decreasing as i = 2,...,d progresses, then b — Yli=i x s l a -s i = 0, and setting 
U := {si, . . . , sj} C a, we have that the columns of Ay are linearly independent. Note that 
U is a cell of A c i and by Lemma l3.7( ii) every maximal cell of A c / arises in this fashion. We 
show that the matrix Ajj has determinant ±1. 

Suppose not. Then the inverse of the matrix Ajj is non-integer, and from the matrix 
equation (Ay)~ l Ay = 1^ we see that there exists a unit vector e,- 6 M. d which is a noninteger 
combination of columns in Ay. Yli=i x s i aL s i = e j- Let z be the vector 

d 

z ■= ^2{x s Ja Sr 

i=i 

Clearly, z € cone(Au) and furthermore z G Z rf since it differs from e, by an integer 
combination of the columns of Ay. So Lemma 13.61 can be applied to b := z: letting 
I := min{i : {x Si } ^ 0} we see that A CjS; (z) = {x Si } < 1 contradicting Lemma 13.61 Hence 
both Ay and (Ay) -1 are integer, their determinant is ±1; since Ay was an arbitrary max- 
imal cell of A c /, we conclude that A c / is unimodular. □ 

The argument concerning the inverse matrix replaces the use of parallelepipeds (compare 
with |35} proof of Theorem 3.1]). Note that all the numbers in the definition of c' are at most 
d d2 , so they have a polynomial number of digits: the perturbed problem has polynomial 
size in terms of the original one, reducing the perfectness test to Problem 13.21 

Remark 3.8. Note that the e of the perturbed vector c' in Theorem 13 . 1 1 made no prescribed 
order on the columns of A. These regular refinements are known as pulling refinements [24] 
and so Theorem 13.11 complements the result of Ohsugi and Hibi [26] regarding pulling 
refinements, of not the normal fan of Q c but, of the polytope Q c itself. See |24] for more 
background on triangulations of polytopes. 

4. Computation 

In this section we wish to provide an illustration of how the results presented in this 
work lead to practical algorithms. We argue that the computational algebra methods for 
detecting the TDI property can be especially efficient in practice when there is a generic 
perturbation of the system in the sense of Proposition 12.91 We wonder if this practical 
efficiency is in some way related to the detection of perfection being in P [9] ? Could 
the perfection recognition problem be solved with a polynomial algorithm based on such 




Figure 1. A chordal graph G with 6 maximal cliques on 10 vertices. 



geometric ideas ? By analogy to the perfectness test, could a unimodular perturbation be 
found in polynomial time (Problem 13. 3p ? 

To show the ideas, we focus on one example of a set packing problem for a chordal graph 
G with 6 maximal cliques on 10 vertices. The coefficient matrix A has 10 rows and 16 
columns and the cost vector c has 16 entries; the first 6 being equal to 1, the last 10 of 
which equal (see the Figure). 

There are essentially five different ways in which we can detect the TDI property in 
the system yA < c coming from the graph G, in finite time, and in polynomial time in 
fix dimension. The first two use Theorem I3.lt Indeed, one can test (even if the time will 
depend exponentially on the dimension) if the triangulation A c / is unimodular. The second 
possibility is to test if the (monomial) toric initial ideal m c i(lA) is generated by squarefree 
monomials, by computing a Grobner basis. 

There are two other ways using Theorem 12.51 Ignoring the generic perturbation arising 
from Theorem 13. 1[ we can simply study the original set packing system yA < c: a third 
algorithm can be based on statement (iii) of Theorem[23J or more precisely, on Corollarv l2.61 
A fourth would be to use statement (v) of Theorem 12.51 

A fifth possibility for the set packing special case is to use the polynomial algorithm 
of Chudnovsky, Cornuejols, Liu, Seymour and Vuskovic [9] which (and actually its early 
predecessors for testing whether a graph is chordal) are certainly the best on this exam- 
ple and any chordal graph, since it provides an immediate positive answer through clique 
separation. However, the example will only contrast the methods expounded upon here. 

The triangulation A c /: Using polymake [15] we can compute the vertices of the simple 
polytope Q c i . We find that there are 288 vertices in all and for each of these vertices, we 
can find their 10 active facets. Finally, for each vertex the coordinates of the active facets 
give rise to a 10 x 10 matrix and we would need to check that each of these 288 square 
matrices have determinant ±1. 

The initial ideal in c '(/y!i): The toric ideal Ia lives in k[a, ...,/, v\, ... , uio] where a, . . . , / 
correspond to the maximal cliques of G (the first 6 columns of A) and where v\ , . . . , v\o 
correspond to the vertices of G (the ordered columns of — iio, the last 10 columns of ^4) as 
before. The toric ideal for set packing matrices A has a very simple generating set: there 
is one generator for every maximal clique in the graph, all of the form zv^v^ ■ ■ ■ v i r — 1, 
where z is the variable corresponding to the clique and «2, • • • , i r } are the vertices of 
that clique. See [271 §2-2] for justification. 
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This simple generating set enables quick computation of the toric initial ideals of A. 
Using Macaulay 2 [19], we computed m c i(lA) = (fv 7 v 8 , ev^VQVio, dv 2 Vg, cv$, bv 3 , av\ ) in 
less than one second on a standard desktop. The monomial toric initial ideal ui c i(Ia) has 
6 generators and all are squarefree. Not only was this initial ideal computed quickly but 
furthermore, because of the small number of generators, we see that the TDI property is 
presented in a highly compact manner, unlike the triangulation A c / above. 

The construction of Ta,c- To construct Ta )C we would first need to find all sets of size 
10 that are not cells in A c . There are a total of = 8008 candidate sets of size 10. 
We would then have to eliminate all 10-sets that are contained in some cell of A c . Again, 
using polymake |15j . we can find all 101 maximal cells of A c . One such cell is the 14-set 
<t := [16]\{8, 13} and so all 10-sets that are contained in a must be removed. Similarly, using 
the other 100 cells of A c we could find all 10-sets that are non-cells. We would then need to 
implement the remaining steps outlined in the algorithm of Corollary 12 .61 a computationally 
challenging process that would need to be carried out for each of the non-cells. This task 
does not compare favourably to the computation of in c /(i^) above. 

The monomial ideal mono(in c /(I^)): There is an algorithm [321 Algorithm 4.4.2] for 
computing the generators of the monomial ideal mono(in c (7 J 4)) which involves four steps, 
all implemented in Macaulay 2 [19]. The first is finding a generating set for the toric 
initial ideal in c (i^); one way to do this is to find a Grobner basis of I a with respect to 
c. For our example, one generating set is {/vyv 8 , ev^v^vio — fv$, dv 2 vg — et>7t>io, cv<± — 
dv^vg, bv 3 — dvQVg, av\ — dv^v^vg} which consists of 1 monomial and 5 binomials in the 
polynomial ring k[o, •••,/, i>i, fro]- Next, we must carry out a multi-homogenization 
procedure on this generating set to form the multi-homogenized ideal in c (I J 4) homo in the 
polynomial ring k[a, ...,/, v\, ... , t>io> A, . . . , F, V\, . . . , V\o] consisting of 32 variables where 
the 16 new homogenizing variables are A, . . . , F, Vi, . . . , Vio- Finding generating sets for 
in c (I^) and in c (/4) homo can be quickly done using Macaulay 2. 

The third step is to find a Grobner basis of in c (/ J 4) homo with respect to an elimina- 
tion order. The fourth and final step would be to extract the monomials in the Grobner 
basis of step 3 - this set of monomials generate the monomial ideal: mono(in c (/ J 4)) = 
{aviv 2 , bcv 2 v 3 V4, bv 2 v 3 v 5 , cv 2 ViVQ, bfv 2 v 3 v 8 , cfv 2 ViV 8 , fv 7 v 8 , dv 2 v 5 v 6 v 9 , dfv 2 v 8 v 9 , aeviv 7 v w , 
bcev 3 v^vjVio, bevsv^vio, cev4VQV?vio, ev§VQV7Vio }. Like the computation of in c '(/^) above, 
this set of 14 generators was done in less than a second on a standard desktop and from it we 
can see the TDI propertry presented in a compact manner. Note that, from the equivalence 
of (iii) and (v) in Theorem 12.51 we also get the 14 k's and r's of Ta,c from these generators. 

We hope that this pedagogical example exhibits some of the utility of computational al- 
gebra methods in integer programming. In larger examples, like Padberg's windmill \29\ last 
figure] (a graph with 21 maximal cliques and 20 vertices), the computation of mono(in c (/ J 4)) 
cannot be carried out on a standard desktop but computing in c i(lA) can be carried out in 
less than a second. However, it should be noted, that for both our example above and for 
Padberg's windmill, it is easy to observe that the graph is perfect via clique separation, a 
classical simple operation that preserves perfectness. 

The ability to only attain in c /(Iyi) in the Padberg example was typical of other examples 
of large set packing problems. Could this or some other "geometric" (polyhedral) method 
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ever handle perfectness - a notion that can be defined in purely polyhedral terms - efficiently 
(in polynomial time) and in a less technical way ? In view of the NP-completeness of the 
general problem [4], [30] , unimodular covering may provide a distinguishing clue for this 
— 1 special case. 
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